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ABSTRACT 

The physical properties of galactic cirrus emission are not well characterized. BOOMERANG is a balloon- 
borne experiment designed to study the cosmic microwave background at high angular resolution in the mil- 
limeter range. The BOOMERANG 245 and 345 GHz channels are sensitive to interstellar signals, in a spectral 
range intermediate between FIR and microwave frequencies. We look for physical characteristics of cirrus 
structures in a region at high galactic latitudes (b~ -40°) where BOOMERANG performed its deepest inte- 
gration, combining the BOOMERANG data with other available datasets at different wavelengths. We have 
detected 8 emission patches in the 345 GHz map, consistent with cirrus dust in the Infrared Astronomical Satel- 
lite maps. The analysis technique we have developed allows to identify the location and the shape of cirrus 
clouds, and to extract the flux from observations with different instruments at different wavelengths and angular 
resolutions. We study the integrated flux emitted from these cirrus clouds using data from Infrared Astronom- 
ical Satellite (IRAS), DIRBE, BOOMERANG and Wilkinson Microwave Anisotropy Probe in the frequency 
range 23-3000 GHz (13 mm 100 /im wavelength). We fit the measured spectral energy distributions with a 
combination of a grey body and a power-law spectra considering two models for the thermal emission. The 
temperature of the thermal dust component varies in the 7 - 20 K range and its emissivity spectral index is in 
the 1 - 5 range. We identified a physical relation between temperature and spectral index as had been proposed 
in previous works. This technique can be proficiently used for the forthcoming Planck and Herschel missions 
data. 

Subject headings: cosmology: observations — cosmology: foregrounds — galactic dust 



1. INTRODUCTION 

Characterising the properties of dust in our Galaxy is an 
important topic of millimeter (mm) and submillimeter astro- 
physical observations. At frequencies above 100 GHz this 
emission is dominated by thermal radiation from large grains 
in equilibrium with the interstellar radiation field. Interstellar 
dust is distributed in filamentary cirrus-like cloud s and cov- 
ers th e sky at both low and high galactic latitudes (lLow et al.l 
[1981 . This emission is usually described by a thermal spec- 
trum, parameterized by the physical temperature of the grains 
Td and by their emissivity versus frequency, which is as- 
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sumed to be a power law with spectral index /?. While 
dust p roperties have been deeply studied i n the galactic 
plane (iDesert et al.l 120081: iDupac et al.l 120031: lLagache et all 
119981) . there is a little i nformation at h i gh galactic lati- 

" 11991 iKiss et al.l 
Dust emission at 



tudes jArendt et al. 1998; Boulange r et al. 
2006; Bot et al. 2009; Pascale et al. 2008)E! 



high galactic latitudes is interesting for two reasons. First, at 
these latitudes the detection of each structure is very little af- 
fected by the overlap of other structures along the line of sight; 
this allows an unequivocal estimation of the physical proper- 
ties of the observed cloud. Second, a good knowledge of dust 
emission at high galactic latitudes is crucial to determine its 
potential contamination of the cosmic microwave background 
(CMB) me asurements, and to improve component separation 
techniques dLeach et al.ll2008l) . 

In the mm range there is a lack of ob servational data, 
and d ifferent models have been proposed by dFinkbeiner et al.l 
I1999L FDS hereafter) to extrap olate the data measured by the 
IRAS dNeugebaue r et al.l 1 19841) to the microwave frequency 
range. 

BOOMERANG -03 offers an unprecedented combination 
of coverage and sensitivity, providing ~ 10° x 10° maps of 
a high latitude region centered at b~ -40°, at 145, 245 and 
345 GHz, with an angular resolution < 10 arcminutes. This 
instrument derives directly from the BOOMERANG payload 
flown in 1997/1998. That payload was recovered and mod- 
ified to make it sensitive to polarization and to improve the 
hardware, keeping almost the same spectral coverage and an- 



19 Knowledge of the dust properties has been improved recentl y by the 
BLAST experiment (see Pascale et al. 2008; Netterfield et al. 2009, and ref- 
erences therein). 
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gular resolution. For furt her inform a tion o n the BOOMER- 
ANG -98 instrument see ICrill et all (120031) and previous re- 
sults on interstellar dust detected by BOOMERANG -98 at 
high galactic latitudes are reported in lMasi etall (1200 1[) . 

In this paper we present an analysis of the characteris- 
tics of diffuse dust emission from far-infrared (FIR) to mi- 
crowave frequencies in the nearby interstellar medium, at 
gal actic latitudes -50 ° < b < -15°, using the BOOMERANG 
03 (iMasi et al.ll2Q06h. Wil ki nson M icrowave Anisotropy Probe 
fWMAP: iHinshaw et al.1 (120071) ). IRAS dNeugebauer et alJ 
119841) a nd Diffuse Infra-R ed Background Experiment 
(DIRBE; iBoggess et al.l (Il992|) ) data. In particular, we focus 
on the observation of eight high-latitude cirrus clouds located 
in the BOOMERANG deep integration field. We derive phys- 
ical parameters of the dust in the clouds, and we study the 
relation between these parameters which can pro vide insight 
into th e nature of the dust grains as suggested by lMeny et al.l 
<2007h . 

The paper is structured as follows. Section [2] describes the 
data sets and the calibration; Section [3] describes the pipeline 
adopted; Section [4] reports the results on the dust properties. 
Conclusions are discussed in Section [5] 

2. DATA PROCESSING 

BOOMERANG 03 (hereafter B03) is a balloon-borne ex- 
periment which in 2003 January performed a 14 day flight 
over Antarctica. It can be considered a pathfinder for the 
High Frequency Instrument (HFI) of the Planck satellite since 
it validated the detectors, the scanning strategy, and initiated 
the relevant data analysis techniques. B03 has observed the 
microwave sky in three frequency bands centered at 145, 245 
and 345 GHz with high angular resolution (~ 10 arcminutes 
at 145 GHz, ~ 7 arcminutes at 245 and 345 GHz). While the 
145 GHz channel is devoted to CMB studies, the two high fre- 
quency channels mainly monitor foreground emission. The 
observed region covers approximately 4% of the sky in the 
southern hemisphere and has been divided in three areas: a 
"deep" (long integration) survey of ~90 deg 2 , a "shallow sur- 
vey" region of ~750 deg 2 and a region of ~300 deg 2 across 
the galactic plane. A detailed description of the B03 instru- 
ment and scanning strategy is published in lMasi etall (|2006|) . 

In this work we study the dust properties by analyzing 
the 245 and 345 GHz channel observations in the deep re- 
gion (70° < R.A. < 95° and -52° < decl. < -39°), which 
provides high signal-to-noise ratio observations of interstel- 
lar dust emission at high galactic latitudes. B03 bands are 
in a particularly key position for the study of interstellar dust 
because they cover the Rayleigh- Jeans part of its spectrum, 
which is currently poorly constrained by observations, and 
they allow us to test the extrapolati on of the models to long 
wavelengths (Finkbeine r et al.ll 19991) . 

In the 245 and 345 GHz detectors, radiation is concen- 
trated on spider web bolometers by cold horns assemblies. 
A metal wire grid is placed in front of each detector, so 
that it is sensitive only to one of the two orthogonal po- 
larizations. Intensity and polarization measurements can 
be derived combining signals from different detectors. In 
this paper we focus on the temperature signal only. Pre- 
vi ous B03 results o n CMB and foreground s are described 
in I Jones et all (120061): IPiacentini et all (120061): iMontrov et all 
(2006); MacTavish et al. (2006); De Troia et al.1 (120071) : 
iNatoli et al.l d2009h : IVeneziani et al.1 ( 120091) . 

In order to study the physical properties of diffuse dust over 



a wide range of frequencies, we include in our analysis the 
IRAS data sets at 100/im (1100 in the following), the DIRBE 
data sets at 240/im (D240 in the following), and the five 5 yr 
WMAP bandit Ka, Q, V and W, centered at 23, 33, 41, 61 
and 93 GHz, respectively. For the IRAS data, we use the Im- 
pr oved Reprocessing of the IRAS Surve y (IRIS £3 described 
in lMiville-Deschenes & Lagachd (120051) . 

2.1. Timeline processing and Cosmic Microwave 
Background removal 

The 245 and 345 GHz channels of B03 (hereafter B245 
and B345, respectively) have four bolometers each, named 
W, X, Y and Z. In this analysis we consider only bolome- 
ters W, Y, and Z at 245 GHz and W, X and Y at 345 GHz - 
the ot her detectors wer e affected by anomalous instrumental 
noise (IMasi et al.ll2006h . 

As in the standard pipeline, the raw time-ordered data 
(TOD) from each detector is deconvolved from its transfer 
function to remove the filtering effects of the readout elec- 
tronics and of the time response of the detectors. The point- 
ing is recovered from on-board attitude sensors, and bad data 
are flagged and not used for the analysis. The deconvolved 
time lines are then reduced with the ROMA map-making 
code (de Gasperis et al. 200 5j) to produce brigh tness maps us- 
ing the HealpbH scheme dGorski et al.l 120051) . Since we are 
looking for large scale structures, we use a pixel size of 13 '.7 
(which corresponds to N S id e = 256) in this analysis. 

At the frequency range and the observed region of the 
BOOMERANG experiment, the dust signal can some- 
times be comparable/subdominant to the measured CMB 
anisotropy. Therefore, to study the physical properties of the 
diffuse dust at the BOOMERANG data, we need first to re- 
move the CMB signal from the 245 and the 345 GHz maps 
before performing any foreground studies on these data. In 
order to remove the CMB signal, we operated as follows: 

1. Since at 145 GHz dust emi ssion is subdomi nant with 
respect to CMB anisotropy (M asi et aLl l2006), we cali- 
brate the BOOMERANG 245 and 345 GHz data using 
the WMAP and the BOOMERANG -98 data. In har- 
monic space, we compute the slope of the linear cor- 
relation between the BOOMERANG 245 (345) GHz 
map and the BOOMERANG -98 150 GHz map, and 
the slope of the linear correlation between the WMAP- 
94 GHz map and the BOOMERANG -98 150 GHz 
map. The ratio of the two slopes is a measurement of 
the calibration factor (in voU/Kcmb) of the 245 (345) 
GHz data of BOOMERANG . For example: 

/ n 245X x 598 \ 
r _ \ a im X a lm I n x 
^ 24 5* " (WMAP ~ n B9S\ W 
\ a tm X a im I 

In other words we calibrate our detectors with respect 
to WMAP, using the B98 map as a transfer to select 
only the CMB part of the map. To account for dif- 
ferent beams and scanning strategy, each a£ m has been 
pr eviously divided b y its window function as discussed 
in IMasi et all (120061) . The calibration factors on CMB 
against the reference channels B245W and B345W are 
reported in second column of Table [T] 
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TABLE 1 

Relative gains of BOOMERanG detectors 



Channel 


1C*w 


1Z*w 


245X 


1.00 ±0.08 


1.09 ±0.02 


245Y 


1.03 ±0.07 


1.15 ±0.02 


245Z 


1.16 ±0.06 


1.17 ±0.02 


345X 


0.85±0.15 


0.92 ±0.03 


345Y 


1.47 ±0.16 


1.26 ±0.02 


345Z 


1.23±0.18 


0.94 ±0.02 



BOOMERANG detectors relative gains for CMB (second column) and dust 
(third column) against the reference channels B245W and B345W, respec- 
tively. The comparison is done with the cross-spectra method (see the text). 
1 - a errors are reported. 

2. We remove the CMB dipole signal from the time lines 
us ing a templa t e deri ved from the parameters given 
in lMather etaLl (11994 . Then we subtract the 145 GHz 
map from the 245 and 345 GHz CMB calibrated time 
lines. 

3. The CMB subtracted time lines are combined into a 
single detecto r map using the ROMA iterative map- 
making code ( de Gasperis et al.l [20051) . The TOD is 
then high-pass filtered at 20 mHz to remove large-scale 
signal from the CMB dipole, system drifts and to cor- 
rect for the non-perfect knowledge of the transfer func- 
tion of detectors. 

4. The relative calibration of the resulting CMB subtracted 
map is corrected for the effect of different spectral 
band-passes by using the IRAS 100 iim map as trans- 
fer: 



245X/245W 



(a 



245X 



x a 



IRAS 

im , 



(a 



245 W v n IRAS\ 
tm Xa £m I 



(2) 



5. These relative calibration corrections have been ap- 
plied to the CMB subtracted time lines to obtain multi- 
detector maps, one at 245 GHz and another at 345 GHz. 
They are reported in third column of Table [T] 

In order to take into account the effect of the BOOMER- 
ANG scanning strategy, the maps derived from the other ex- 
periments are projected onto four time- streams correspond- 
ing to the four BOOMERANG photometers. The time lines 
are then filtered in the same way and assembled in a multi- 
detector map replicating the flight pointing and flagging of 
BOOMERANG . We correct for the slight effect of this pro- 
cedure on the measured fluxes as described in Section 3.1. 
The CMB from the WMAP maps has been removed using 
the 145 GHz BOOMERANG channel. All our maps are 
thus CMB cleaned, with the assumption that the BOOMER- 
A NG 145 GHz data are CMB dominated; indeed, in Fig. 33 
of iMasi et all (120061) it is shown that at 145 GHz dust rms is 
much less than CMB rms. Nevertheless, it is possible that 
some residual dust flux is present in our 145 GHz map in the 
cirrus clouds regions. We estimated the possible dust resid- 
ual using the recovered models (see Section 4) and we find 
that the contamination can be of the order of 5% of the signal 
at 345 GHz expressed in CMB temperature units. This con- 
tamination on brightness has the same value at all frequencies 
when expressed in CMB thermodynamic units, and we have 
seen that it amounts to a fraction of the error bars of the mea- 
sured average brightness of the cirrus clouds at all frequen- 
cies. 



3. DATA ANALYSIS 

The goal of this work is the estimation of the physical prop- 
erties of cirrus structures through the measurement of their 
flux over a wide frequency range. Thus, when measuring the 
flux, we need to collect all and only pixels with dust emission. 
Since dust emission peaks at high frequency, we choose the 
maps with the highest signal-to-noise ratio, B345 and 1100, 
to select the brightest pixels belonging to dust structures. We 
then combined B345 and 1100 in pixel space producing a cor- 
relation map given by: 



M{p)i 



IB345(P) X Inoo(p) 



N2 J2p' ^5345 0?') ^Z p ' I'i 
IB345(P) X IIIOO(P) 



noo(P') 



^5345 X CT/100 



(3) 



where p identifies a pixel in the Healpix scheme, / is the in- 
tensity map, and a is the rms of the intensity in the observed 
region. The equality in Equation [3] holds since the average 
of the intensity maps is vanishing because of the high-pass 
filtering procedure. The normalization has been chosen in 
such a way that the correlation between two identical maps, 
summed over all pixels, is 1. Since we are interested in posi- 
tive emission regions only, we masked the negative emission 
ones when computing the correlation map. 

With this procedure we could identify eight large clouds 
which are marked by black circles in Figure [T] As a starting 
point for the pixel selection, in each of these areas we select 
the brightest pixels by imposing a threshold M(p) > 1 .5. This 
value allows us to include cold sources that have a faint signal 
in the 1100 map. In order to preserve the irregular shape of 
the clouds, later on we increased the initial area in steps, ex- 
tending it by 0.1° around each pixel. We choose this step size 
since pixels have ° .22 side (N s id e = 256), so 0°.l starting 
from the center of the pixel is the minimum distance to col- 
lect another pixel. We patch a maximum distance of 1° from 
the starting mask, to avoid contributions from other structures. 
In this way, we preserve the shape of the significantly cor- 
related region in the correlation map, and include the other 
pixels dominated by dust emission, without missing any. The 
integrated flux increases and reaches a maximum; then it de- 
creases if we further widen the selected area. In fact, the de- 
tectors are AC coupled, and produce negative bounces after 
detecting a positive signal from a localized source. The error 
bars on the flux are given by adding in quadrature the errors of 
each pixel belonging to the mask. We use this effect to detect 
the area to be included in the flux measurement deriving the 
dimension of the mask corresponding to the maximum flux 
through a Gaussian fit. We then associate with each cloud the 
flux and the error corresponding to that mask size. The bright- 
ness maps are shown in Figure [2 and flux values and errors 
are reported in Tabled An example of this technique is shown 
in Figure [3] and the increasing mask method is shown in Fig- 
ure 01 From Figure [3] we can see that the flux does not depend 
strongly on the mask radius. In WMAP Q, Ka and K bands 
the beam size is comparable to the size of the smaller cirrus 
clouds. In these cases we measure the flux as the integral of 
the brightness over a circle centered in the source coordinates, 
with a 3ctbeam radius. We also check that this flux is consis- 
tent with the flux measured with the method described above. 

We apply this procedure for each cirrus cloud at each fre- 
quency taking into account the noise and beam characteris- 
tics of each experiment. Pixel standard deviation has been 
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TABLE 2 

Clouds locations and measured flux 



Region # 


i 
i 




Q 
J 


4 


c 

J 





1 


Q 

o 


Ra (deg) 




o4. / 


OO A 

50.4 


oU.j 


no o 


O /.U 


oU.U 


78.5 


O/l c 

o4.j 


Dec (deg) 




zLS ^ 


-'H-o. 1 


4S S 


-'HO.U 


40 S 


4Q n 


49 n 


4Q S 


Area (deg 


2 ) 


0.4 


5.8 


8.0 


2.7 


5.6 


6.0 


2.5 


0.5 


s K [Jy] 




<0.1 


<0.1 


3.7 ±0.4 


<0.1 


0.7 ±0.1 


<0.2 


<0.2 


<0.2 


s Ka [Jy] 




<0.2 


<0.2 


2.1 ±0.2 


<0.1 


0.8 ±0.2 


<0.2 


<0.3 


<0.2 


Sq 


Jy 






<0.3 


<0.4 


1.1 ±0.2 


0.2 ±0.1 


0.8 ±0.2 


<0.3 


<0.4 


<0.6 


Sv 


Jy 






<0.5 


<0.6 


<0.5 


<0.7 


1.7 ±0.7 


<0.9 


<0.7 


<0.7 


Sw 


[Jy 






< 1.0 


< 1.0 


< 1.0 


< 1.1 


1.7 ±1.3 


< 1.3 


< 1.0 


< 1.1 


SB245 


Jy 




<2 


14±5 


8±4 


3±2 


7±6 


26 ±12 


6±2 


0.6 ±0.4 


Sb345 


Jy 




5±1 


83 ±33 


90 ±33 


7±2 


74 ±20 


90 ±32 


22±8 


3.5 ±1.6 


SD240 


Jy 




81 ±30 


778 ±385 


500 ±354 


384 ±175 


686 ±224 


430 ±345 


187±81 


72 ±29 


Shoo [Jyl 




44.2 ±5.3 


560 ±102 


585 ±151 


185 ±25 


446 ±86 


305 ±95 


99 ±16 


<8 


N(HI) 


[I0 21 cm- 2 1 


0.375 


0.622 


0.374 


0.433 


0.441 


0.222 


0.323 


0.331 



Coordinates of the centers and areas of circled regions in Figures [T] and [2 (second, third and fourth row, respectively). The following ro ws r eport flux values 
registered by 5 yr WMAP, BOOMERANG , D240 and 1100 channels corrected from the bias with the procedure described in Section [3~T1 In the last row 
corresponding HI column densities are reported. 



estimated from noise maps for BOOMERANG and WMAP. 
BOOMERANG standard deviation also takes into account 
also calibration errors, which increase the flux error bars by 
~ 10%. In the case of IRAS and DIRBE we computed the rms 
of the data in several "clean" regions at high galactic latitudes 
and took the average standard deviation of the data as pixel 
error. 

In Figur^H the areas that we have identified as cirrus clouds 
are circled, and the approximate coordinates of the centers, 
as well as the cirrus areas, fluxes and column densities, are 
reported in Table [2 The hydrogen column densities have 
bee n measured using the Leiden/ Argent i ne/Bo n n (LAB) sur- 
vey iKalberla et alJ (120051) : lArnal et alJ d2000h : iBajaja et alJ 
d2005h . 

We compared the BO OMERANG fluxes to th ose predicted 
by the model 8 of FDS dFinkbeiner et al] 1 19991) at the nom- 
inal BOOMERANG frequencies, measuring the fluxes from 
extrapolated simulations of dust maps in the clouds regions 
we identified. A flux-to-flux comparison, made by plotting 
the real data fluxes versus the predicted ones, has a best fit 
line with slopes 1.5 ± 0.3 and 0.9 ± 0.3 at 345 and 245 GHz 
respectively. 

3.1. Estimation of cirrus physical parameters 

It is generally accepted that the submillimeter spectrum 
of thermal dust can be expressed as a blackbody times a 
frequency dependent (power-law) emissivity factor. Simple 
emission models predict the emissivity factor to be approxi- 
mately 2. Significant variations to the (3 = 2 emissivity law 
occur when taking into account the diso rdered structure of 
amorphous dust grains (Meny et al. 2007). At tens of arcmin- 
utes angular scales, it is also possible to detect different tem- 
peratures in the structures due to different dust com ponents as 
for example in model 8 of Finkbein er et al.l ([1999). In these 
cases, an isothermal model could lead to spectral index values 
smaller than 2 while a two temperatures model with a cold 
dust component at approximately 10 K and (3 = 2 could better 
approximate the physics of the emission. Both these models 
are tested on our sample. 

In order to analyze cirrus characteristics (e. g., tempera- 
ture and emissivity) we fit the flux data taking into account 



two contributions: a blackbody times a frequency-dependent 
emissivity emission component, peaking in the far infra-red 
region, and a second emission component at lower frequen- 
cies parameterized as a power law with a spectral index able 
to reproduce synchrotron or free-free emission. When assum- 
ing the thermal dust emission to be isothermal, the spectral 
energy distribution (SED) results in 

S(v) = A 1 (^J +A 2^) bb(^,T d ) (4) 

where bb(^,T d ) is the standard blackbody function. u and 
v\ are 33 GHz and 3000 GHz, respectively. Td is the tem- 
perature of the considered dust cloud, (3 is the corresponding 
spectral index, and a is the spectral index of the low frequen- 
cies contribution. A\ and A 2 represent the amplitude of the 
two components, which at these galactic latitudes is propor- 
tional to the optical depth. When modeling the thermal dust 
emission with a warm and a cold component we set the spec- 
tral index of both to the theoretical value of 2 and the SED 
becomes 

(5) 

where T c , T w are the temperatures of the cold and warm com- 
ponents respectively and B 2 and B 3 the corresponding ampli- 
tudes. T c is set to 10 K while T w varies. Bi and a 2 describe 
the amplitude and the emissivity of the low frequency compo- 
nent, respectively. 

From the measured SED we estimate the parameters A\ , A 2 , 
a, (3 and Td in the isothermal model, and the parameters B\, 
B 2 , #3, T w , a 2 in the two components mod el, using a Monte 
Carlo Markov Chain (MCM C) algorithm (IChristensen et al.l 
1200 ll: Eewis & Bridle! 12002). Following Bayesian statistics, 
the probability to have a set of parameters p given a set of 
data d is 

P(p|d)cxP(p)P(d|p), (6) 

where P(p) is the a priori probability density of the parameters 
set p and ,P(d|p) is the probability density of data d given a 
set of parameters p, named the likelihood function. The poste- 
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TABLE 3 
Priors 



Correlation map B345 x 1100 
0.00c 




80 82 84 
RA [Deg] 



Fig. 1. — Correlation map (Equation[3} between B345 and 1100. All posi- 
tive structures shown in the map above were included in the final mask, with 
the exception of the regions centered at (R.A., decl.) = (78°, -40°), and (88°, 
-50°) and the region 85° < R.A. < 90° and -46° < decl. < -42°. They are 
small parts of very complex shaped structures which have their largest part in 
the BOOMERANG shallow field. 

rior probability P(p|d) is estimated using MCMC algorithm. 
Given a set p;, with likelihood Li and posterior probability 
P(p-|d), the MCMC algorithm generates an independent set 
Pz+i with likelihood L i+ \ and posterior probability P(p /+1 |d). 
This second set is accepted according to a rule which also 
guarantees a good sampling of the probability density in a rea- 
sonable computational time. For example in the Metropolis- 
Hastings algorithm the new set p/+i is always accepted if 



A(7+l,0 = 



P( P/+1 |d) L; +1 P( P . +1 ) 



> 1 



(7) 



This guarantees the convergence to the maximum of the likeli- 
hood function. If instead A(/+ 1 , i) < 1 , the new set is accepted 
with a probability proportional to the ratio A(/+ 1 , /). This en- 
sures a good sampling of the distribution. F or a better descrip- 
tion o f the operation of the algori thm see IChristensen et al.l 
(120011) and iLewis & Bridlel (120021) . We chose wide flat a 
priori probability densities of the parameters, as reported in 
Tab. [3l Three chains for each cloud have been run, starting 



Parameter 


Limits 


Oi. 


free 


T d 


0-40 


T w 


0-40 


P 


-1-10.0 


Ai 


>0 


A 2 


>0 


a 2 


free 


Bi 


>0 


B 2 


>0 


B 3 


>0 



List of a priori probability density imposed on parameter set. 

from different points of the parameter space. The convergence 
has been checked using the test algorithm encoded in the Get- 
Dist software, available in the public CosmoMC packagS 

The combination of the BOOMERANG scanning strategy 
and AC coupling induces an effective filtering to the sky sig- 
nal. To take the induced bias into account, we applied the 
same effective filtering to all the other data sets as described 
in Section [2 and estimated its effect by means of simulations. 
Since filtering effects depend strongly on the scanning strat- 
egy, which translates in the map into a dependency on the 
position of the cloud and on the nearby structures, we pro- 
duced for each cirrus cloud 100 simulations with the same 
shape, position and noise properties of the measured cloud, 
and processed each simulation through the pipeline applied 
to the data. Comparing the average output flux values with 
the input ones we obtain the bias factors the measured fluxes 
have to be corrected for. The factors show that approximately 
20% of the sky signal is dumped by the pipeline. As expected, 
these factors depend on the dilution of the cirrus in the beam 
and on position and shape of the cloud itself. 

4. RESULTS 

Our analysis shows that the two temperatures model of 
Equation [5] and the isothermal model of Equation [4] are sta- 
tistically equivalent. However, the 10 K component is poorly 
constrained in almost all clouds, as its amplitude B3 is con- 
sistent with zero. The temperature of the warm component 
T w is consistent with the temperature Td estimated with the 
isothermal fit. This consistency provides a further check on 
the temperature values obtained. 

We can distinguish two sets of cirrus clouds: the first con- 
sists of two of the eight selected dust regions (marked with 
numbers 3 and 5 in Figure [U whose emission is clearly de- 
tected in the full range of analyzed frequencies; we are then 
able to derive both the thermal dust and low frequency com- 
ponent. The SED for these clouds, estimated with both the 
isothermal and the two-temperatures model, is shown in Fig- 
ure The remaining six clouds (numbers 1, 2, 4, 6, 7 and 8) 
form the second set and present the thermal dust component 
only, as plotted in Figured This can be seen in the values of 
the parameter A\, and consequently of the spectral index a, 
which in the second set are un-determined. 

The characteristics of each cloud for the isothermal and 
the two temperatures model are reported in Tables |4] and 
respectively. In order to check the goodness of the fit, the 
last columns of these tables report the cumulative distribution 

23 |http : //c osmologist . inf o /cosmomc/1 
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FIG. 2.— The BOOMERANG deep region observed, from top bottom and from left to right, by WMAP at 23, 33, 41, 61 and 93 GHz, by BOOMERANG at 
245 GHz and 345 GHz, by DIRBE at 240 /xm and by IRAS at 100 /im. 
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0.0 0.2 0.4 0.6 0.8 

Mask size (cleg) 



1.0 



Fig. 3 . — Plot of the integrated flux vs. size of the integration area for cirrus 
2 at 345 GHz. The initial region, corresponding to 0°.0, is selected using the 
correlation map (Figure 1) and then increased with the technique described 
in the text in order to preserve the shape of the cloud. The solid line is the 
best Gaussian fit on the estimated fluxes. In this case, the maximum flux is 
found in a mask increased of 0° .45 with respect to the starting one. The final 
value of the flux and associated error are those measured with this mask size. 

Mask size = 0.0/0.25/0.45 deg 
0.00 

0.0 0.1 0.2 0.3 0.4 



-48 7 




88 
RA [Deg] 

FIG. 4. — Starting mask (inner contour) and 0° .25 and 0° .45 "radius" mask 
(medium and outer contour respectively), the latter corresponding to maxi- 
mum integrated flux, for cirrus 2 at 345 GHz (same as in Figure 3). 



function P of a x 2 distribution with the proper number of de- 
grees of freedom. The ideal fit should have P ~ 0.5; P « 0.5 
is an indication of overestimated error bars; P » 0.5 indi- 
cates that the fitting function is not a good model. 

One and two-dimensional posterior probability of parame- 
ters in Equation|4] are shown for the cloud number 3, as exam- 
ple, in Figure [7] 



Assuming an isothermal model of dust emission the clouds 
have dust temperatures in the range 7^ G [6; 21] K, and emis- 
sivities /3 G [1.1; 5]. The amplitude and spectral index of the 
warm dust and the low frequency component estimated us- 
ing the two components model are fully consistent with the 
correspondent values estimated with the isothermal model. 

In order to study dust emissivity variation across the patch 
and its dependence on dust temperature, in the following anal- 
ysis we will refer only to the isothermal model for the dust 
emission. 

4. 1 . Temperature dependence of the spectral index 

Dust emissivities and temperatures in the BOOMERANG 
deep field cover a range of values (see Table 0]). The spectral 
shape of the thermal component of the dust emissivity (second 
term of Equation [4]) is strongly dependent on a combination 
of and /?, whose two-dimensional posterior probability has 
a elongated, slant shape (Figure Q, indicating a degeneracy 
between these two parameters. Two more degeneracy figures 
are obtained in the A 2 -f3 and A\-a planes, as e xpected from 
the fu nctional form. Previous works on P ronaos (iDupac et alJ 
120031) and Archeops dDesert et al.l 120081) data which analyze 
many cold clouds find an inverse relation between the dust 
temperature and its spectral index, the former being expressed 
as a function of the latter. 

Given the shape of the posterior probability we find in our 
data, we want to investigate whether in our case the /3-T& re- 
lation is a physical characteristic of the dust or just a conse- 
quence of the functional form and of measurements errors. 
We run then a MCMC on each region expressing, in Equa- 
tion HI f3 as a function of Td following the two models: 



1. iDesert et al.1 (120081) . defined by the relation 

P=AxT d p 



2. IDupac et all (120031) . through the relation 

1 



(8) 



(9) 



where A = 1 1.5 ± 3.8, p = -0.66 ± 0.05, S = 0.40 ± 0.02, uj = 
0.0079 ±0.0005 K~ l are the values derived from the Archeops 
and Pronaos data, respectively. 

In order to check the concordance between B OOMERanG 
and the other considered experiments, we fit our sample us- 
ing the two models quoted before (Equations 8 and 9) includ- 
ing the error induced by the degeneracy shape relation. We 
thus perform the fit on eight random points, one for each cir- 
rus, within the 68% contour of the two-dimensional posterior 
probabilities. We repeat this procedure 60000 times for each 
model, a number of steps which allows a good sampling of the 
distributions of the parameters A, p, 5 and uj in Equations 8 
and 9 respectively. This technique allows a complete sam- 
pling of the degeneracy, not being dominated by its shape in 
the final result. From the peaks and the standard deviations of 
the distributions we derive as best fits: 



P = (33.3 ± 6.0) xT d - 11±01 



P 



1 



model 1 
model 2 



(0.035 ± 0.004) x T d - (0.013 ± 0.030) 

finding a trend similar to the Archeops result. Figure [8] shows 
posterior probabilities 68% contour plot in the Td-/? plane, 
for all the analyzed structures together. The dashed and solid 
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TABLE 4 

Physical characteristics for each dust structure in the iso-thermal model 



Region # 


a 


T d [K] 


P 


logAi 


logA 2 


P(S) 


1 




17.0 ±2.3 


1.8±0.3 


<-1.5 


-2.2 ±0.4 


0.52 


2 




16.9±2.8 


1.8 ±0.4 


<0.1 


-1.2±0.5 


0.79 


3 


-2.0 ±0.2 


20.4 ±5.4 


1.5±0.5 


0.3±0.1 


-1.9±0.6 


0.71 


4 




15.1 ±1.2 


2.6 ±0.3 


< 1.6 


-1.2±0.3 


0.53 


5 


0.4 ±0.3 


19.0±2.8 


1.4±0.3 


-1.4±0.2 


-1.7 ±0.4 


0.69 


6 




18.0±3.9 


1.3±0.3 


<-1.8 


-1.6±0.5 


0.91 


7 




19.6±4.1 


1.1 ±0.5 


<-0.2 


-2.4 ±0.5 


0.47 


8 




6.5 ±2.6 


5.1 ±1.8 


<-3.5 


1.7±1.3 


0.56 



Physical parameters of the clouds assuming an iso-thermal (Equation[4j model of dust emission, as measured using an MCMC algorithm. The errors correspond 
to 68% confidence interval. 



TABLE 5 

Physical characteristics for each dust structure in the two component model 



Region # 


OL 2 


T w [K] 


logBj 


logB 2 


logB 3 


P(S 2 ) 


1 




17.1 ±2.6 


<-4 


-2.3 ±0.4 


<-2 


0.62 


2 




18.7±4.7 


<0.1 


-1.6±0.6 


<-0.8 


0.82 


3 


-2.2 ±0.2 


21.8 ±5.4 


0.3±0.1 


-2.0 ±0.6 


<-0.9 


0.65 


4 




17.7 ±0.5 


0.3±0.1 


-1.7±0.1 


<-4.5 


0.56 


5 


0.4 ±0.2 


20.3 ±5.5 


-0.2 ±0.1 


-1.8 ±0.7 


<-0.6 


0.68 


6 




17.6±4.5 


<-0.4 


-1.7±0.8 


-0.4 ±0.3 


0.76 


7 




21.1 ±3.1 


<-0.2 


-2.7 ±0.7 


-1.0 ±0.4 


0.38 


8 




9.0±2.1 


<-2.3 


-1.7 ±0.7 


<-0.1 


0.62 



Physical parameters of the clouds assuming a two components (Eq.[5} model of dust emission. The errors correspond to 68% confidence interval. The amplitude 
and spectral index of the low frequency component and of the warm dust component are fully consistent with the ones estimated using the isothermal model. The 
cold component, described by the coefficient B3 is determined in two clouds. 



Region 3 Region 5 




0.1 1 , , , I 0.1 1 , , , 

10 100 1000 10000 10 100 1000 10000 

Freq (GHz) Freq (GHz) 

Fig. 5. — Spectra of four cirrus clouds from 23 to 3000 GHz (13 mm to 100 fim in wavelength). Solid line: fitting curve of Equation|4] Dashed line: fitting 
curve of Equation \5\ Low frequency emission is detectable and is fitted by means of a frequency dependent power law. 



lines show BOOMERANG parameters using model 1 and 2 
respectively. The dotted and dash-dotted lines mark the best 
fits obtained by Archeops and Pronaos experiments respec- 
tively. 

4.2. Emission at low frequencies 

As previously discussed, a low frequency emission is de- 
tected in regions 3 and 5. The spectral index a obtained from 
the fit of the flux in Jy can be converted in the correspond- 
ing spectral index 7 in antenna temperature by the relation 

7 = a — 2. 

a is the spectral index reported in the second column of Ta- 
ble |U According to this conversion we obtain 7 = -4.2 ± 0.2 
in region 3 and 7 = -1.6±0.2 in region 5. We expect the 
detected emission to be a combination of different emissions 



dominating in these bands, typically free-free and synchrotron 
with their typical power laws. Since the considered fre- 
quency range is limited we approximate the combination of 
the two emissions with a single power law. We also esti- 
mate the contribution of these components including the Ha 
SH ASSA (The Southe rn H- Alpha Sky Survey Atlas £3 sur- 
vey dFinkbeinerl |2003|) as mo nitor of free-free em ission and 
the Haslam 408 MHz survey (lHaslam et al.Hl 9 82h as monitor 
of synchrotron emission. Estimated fluxes are reported in Ta- 
ble [6l these values have to be considered as upper limits on 
the emission. Galactic free free is negligible in both regions 
while synchrotron is dominant. 



http ://amundsen. swarthmore.edu/ 
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Region 4 Region 6 
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Fig. 6. — Spectral energy distributions of three cirrus clouds from 245 to 3000 GHz (1.2 mm to 100 /xm in wavelength). Solid line: fitting curve of Eq. |4] 
Dashed line: fitting curve of Eq.[5] These sources don't have a significant low frequency emission so only the thermal dust emission is shown. 



TABLE 6 

Free-free and synchrotron contributions 



Region # 


K 


Ka 


S ff [Jy] 
Q 


V 


w 


K 


Ka 


S s [Jy] 
Q 


V 


w 


3 
5 


<0.1 
<0.2 


<0.1 
<0.1 


<0.1 
<0.1 


<0.1 
<0.1 


<0.1 
<0.1 


<2.1 
<0.3 


<0.9 
< 0.4 


<0.4 
<0.2 


<0.1 
<0.1 


<0.1 
<0.1 



Upper limits on free-free and synchrotron emission in WMAP bands, estimated in regions that show a non negligible flux at low frequency. 



4.3. Spinning dust 

There has been r ecent interest in spinning dust emission 
both theoret ically (Draine & Lazarian 1998) and observa - 
tionally (e.g. lde Oliveira-Costa et al.ll2008l: iGold et alJl2009h . 
Large dust grains in cirrus clouds have a thermal emission 
which peaks at ~ 250 /xm, but smaller grains emit transiently 
at higher frequencies, and could also produce rotational emis- 
sion by spinning at microwave frequencies. This could be 
relevant below 60 GHz. 

The only regions showing a non-negligible low frequency 
emission are 3 and 5. Region 3 includes a radio source that 
dominates the signal at low frequency so we do not expect to 
find a significant spinning dust detection. In order to check 



whether signal in region 5 can be generated by rotational dust 
we have to remove free-free and synchrotron contributions, 
reported in Tabled from the measured fluxes before estimat- 
ing the spinning dust level. After the removal, no detection is 
found in that cloud either. 

5. CONCLUSIONS 

This paper analyzes the properties of galactic cirrus clouds 
at high latitude combining the high-frequency BOOMER- 
ANG dataset together with IRAS, DIRBE and WMAP. We 
located eight clouds in the deep survey area of BOOMER- 
ANG and for each of them we estimated dust temperature and 
emissivity spectral index finding the temperature to be in the 
7 - 20 K range and its emissivity spectral index in the 1 - 
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Fig. 7. — One and two-dimensional posterior probability of the parameters for one of the clouds. The posterior probability is obtained using an MCMC 
algorithm. The black ellipses show the 68% and 95% confidence levels. In this case we consider region 3, and parameters values are a = -2.0 ±0.2, Td = 
20.4 ±5.4, p = 1.5 ± 0.5, log(Ai) = 0.3 ± 0.1, log(A 2 ) = -1.9 ± 0.6. 




Fig. 8. — Two-dimensional 7j-/3 posterior probability contour plot of the 
whole structures set. The contours are at 68% of confidence level; the points 
inside the contours are used for fitting Equations 8 and 9 as explained in 
the text. The solid an d da shed lines show BO OMERANG best fit follow- 
ing Desert et al. ( 2008) and Dupac et al. (2003) models respectively. The dot 
and dot-dashed lines show the best fits from Archeops and Pronaos datasets, 
respectively. 



5 range. We also investigated the possibility of a two tem- 
peratures model to describe the dust emission, having both 
the components (3 = 2 and the cold component temperature 
T c = 10 K. Within this model the temperature of the warm 
component is consistent with the temperature estimated using 
the isothermal model within la error bars while the cold com- 
ponent is detected in regions 6 and 7. This, and the presence 
of a 6.5 K region, shows the presence of cold dust at high 



latitudes. 

Taking into consideration the shape of the joint poste- 
rior probability of the parameters estimated with the isother- 
mal model for each observed object, our data confirm a 
model in which temperature and spectral index are in- 
versely corr elated, as was sug gested in previ ous analyses 
on Pr onaos dDupac et al] 120031) and Archeops dDesert et al] 
[2008|)data. 

A comparison to the e xtrapolation of IRAS data 
from Finkbe iner et al] (119991) indicates that this predic- 
tion has limited accuracy in the range of frequencies observed 
with BOOMERANG . In particular, the model seems to 
underestimate the dust brightness at 345 GHz by a factor of 
~ 1.5, being the equality with the model within 2a. 

At lower frequencies we detect a signal in the observed re- 
gions which is mostly generated by synchrotron and free-free. 
We do not find any evidence of spinning dust. 

From this analysis it is clear that new data are required to 
improve the knowledge of properties of high galactic latitude 
cirrus clouds, and of dust in general. The method developed 
in this paper allows to identify the location and shape of dust 
clouds, and to extract the flux from observations with differ- 
ent instruments at different wavelengths and angular resolu- 
tions. This technique can be proficiently used to analyze the 
forthcoming Planck and Herschel data sets, which will pro- 
vide higher sensitivity, wider spectral range and, in the case 
of Planck, full sky coverage. 
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